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Tree-level and complete one-loop parametrisation of the linear sigma model (LcxM) is 
performed and the phase boundary between first order and crossover transition regions 

O 



of the m w — m^-plane is determined using the optimised perturbation theory (OPT) as 
a resummation tool of perturbative series. Away from the physical point the parameters 
of the model were determined by making use of chiral perturbation theory (ChPT). The 
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^3 ■ = were estimated. 



location of the phase boundary for m w = and of the tricritical point (TCP) on the 



^ . 1. Introduction 

sh ; 

Lattice studies show that at the physical point of the m u ^ — m s -plane the transition as 
a function of the temperature is of crossover-type [[I]. In the m u ^ = rn s = limit effective 
models predict a first order phase transition and for m s = oo and m u ^ = a second order 
phase transition [|2]. Theoretically it is interesting to know where the boundary of the 
first order region lies and where the TCP is located on the m u = axis. 

Several lattice studies with degenerate quarks (m u ^ = m s ) show that the critical pion 
mass on the boundary decreases substantially, from m£ ~ 290 MeV to m% = 67(18) MeV, 
as finer lattices and improved actions are used (see e.g. [EJI1] ). A low pion mass allows for 
studies in effective models based on the chiral symmetry because it is expected that they 
work better closer to the chiral limit. Existing works performed in these models give also 
a low value for the critical pion mass [ El El El • in these works only the Goldstone masses 
were tuned, the other parameters of the Lagrangian were kept at their values determined 
at the physical point. We think that in this low mass regime effective model studies might 
complement lattice investigations if coupling parameters can be determined accurately. 
Hence we let the coupling constants to vary with the pion and kaon masses when moving 
away from the physical point and determined their values using the low energy results of 
the ChPT, namely the dependence of f n , fx,^ and on m n and tuk- 
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2. The parametrisation of the linear sigma model 

The Lagrangian is written in terms of a 3 x 3 matrix M = A a (<r a + «7r a )/v / 2 (A a : a = 
1 . . . 8 are the Gell-Mann matrices and Ao := a/2/31) containing the nonet of scalar (er a ) 
and pseudoscalar (n a ) particles: 

L(M) = ^Tr (d^d^M + li 2 M^M) - f x (Tr (M ] M)f - f 2 Tr (M f M) 2 

-g (det(M) + det(M t )) + e a + e 8 a 8 . (1) 

The determinant term breaks the axial U(l) symmetry and the external fields are in- 
troduced in order to give masses to the pseudo-Goldstone bosons. From the original 
0-8 basis one can switch with an orthogonal transformation to a basis in which in the 
broken symmetry phase the two non-vanishing expectation values are the non-strange 
(x = (v2(cro) + (o"s))/\/3) and the strange (y = (((To) — y/2(a&)) / condensates. 

2.1. Tree-level parametrisation 

The tree-level parametrisation of the model can be performed at zero temperature by 
solving a set of coupled linear equations [ |HJ E] • Only the trace of the mass matrix is 
used in the mixing r\ — rf sector. Unfortunately the values the quartic coupling fi and 
the mass parameter ii 2 can be obtained using only the scalar sector which compared to 
the pseudoscalar sector is experimentally less known. Besides, one knows nothing about 
the and m^-dependence of the scalar masses and as a consequence some assumptions 
are needed. For example in [ |5] we required the Gell-Mann-Okubo mass formula to be 
satisfied in the scalar sector on the entire m n — m^-plane with the same accuracy as at 
the physical point. Performing a consistency check by comparing m v and m v ' calculated 
in the LaM with their values given by the ChPT remarkable agreement was obtained up 
to a limiting value of mx ~ 800 MeV, even away from the physical point. 

2.2. One-loop parametrisation 

The parametrisation at the one-loop level of the perturbation theory is more compli- 
cated and there are many possible ways of selecting the set of non-linear equations. One 
such set of equations is presented below (see [[TU] for details). 

Because in the broken symmetry phase the tree-level mass can be negative a resumma- 
tion is needed. We performed this using the optimised perturbation theory method of Ref. 
[ E] where a mass term is added to and subtracted from the Lagrangian, the difference 
between the original and the new mass is treated perturbatively (first at one-loop level). 
The new mass is determined from the principle of minimal sensitivity which requires that 
the one-loop pole mass for the pion (M v ) stays equal to its tree- level value (m^). This 
can be translated into a self-consistent equation for the pion mass 

ml = -ix 2 + (4/ x + 2/ 2 )x 2 + 4f lV 2 + 2gy + ReZ n (p 2 = m 2 n ). (2) 

Two more equations are the one-loop pole mass of the kaon and rj (smaller mass eigenvalue 
in the mixing sector): 

M 2 K = -Li 2 + 2(2/ x + f 2 )(x 2 + y 2 ) + 2f 2 y 2 - V2x(2f 2 y - g) + ReL K {p 2 = M 2 K ), (3) 
(p 2 - ml x - E Vxx (p 2 ) -m 2 xy - E Vxy (p 2 ) \ 
V - ^ y (p 2 ) P 2 - ml y - H Vyy {p 2 ) ) p2=M ^ ~ °- W 
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We require additionally the equality of the one-loop an tree-level kaon masses: 

M 2 K = m\ = ml- 2gy + Af 2 y 2 - V2x(2f 2 y - g), (5) 

and also make use of the PCAC relations: 

U = Z~^M-\-iD-\ V = 0))x, f K = Z K ^M K 2 {-iD K \p = 0))(x+v / 2y)/2. (6) 

Away from the physical point we use for the determination of m v , f n and fx the formulae 
of the SU(3) ChPT in the large-iV c limit. Because in the one-loop parametrisation the 
renormalisation scale (7) appears as a new parameter we checked how the other parameters 
and masses depends on it. The wave function renormalisation constants have a plateau 
and the variation of the one-loop masses is the mildest in the range / e (1000, 1400) MeV, 
which was selected to be used for the determination of the phase boundary. 

3. Thermodynamics of the LctM 

The order of the phase transition is determined by solving simultaneously the self- 
consistent gap equation for the pion mass 

ml = -u 2 + (4/x + 2f 2 )x 2 + Af lV 2 + 2gy + £> 2 = m 2 , m^m,), I), (7) 

where T ln (p,m i (m 7T ),l) is the pion self-energy, and the two equations of state 

e x = -fx 2 x + 2gxy + Af 1 xy 2 + (Af 1 + 2f 2 )x 3 + j: i J t tf(x,y)l tp [m l (m 7T ),T}, 

e y = -fx 2 y + 2gx 2 + Af 1 x 2 y + (Af 1 + Af 2 )y 3 + j: i J l t y t (x,y)l tp [m t (m w ),T]. (8) 

The isospin multiplicity factor Jj and the coefficients tf and if are given in Ref. [Ej. 

The difference between the tree-level and one-loop parametrisation is that in the first 
case one are enforced to throw out entirely the vacuum fluctuations and use only the 
temperature dependent part of the integrals, while in the second case the renormalised 
vacuum fluctuations are taken into account. The difference between the two methods of 
solving the model was checked quantitatively in [llOj. 

4. The phase-boundary in the m n — mx— plane 

In the case of the tree-level parametrisation the result on the phase boundary is not too 
conclusive because of the dependence of the phase boundary on the assumptions made 
for the scalar sector in the process of parametrisation. We have not seen the trace of 
the tricritical point. We estimated that the phase boundary crosses the diagonal of the 
m w — m^-plane for critical pion mass in the range m^. = 40 ± 20 MeV. 

As for the one-loop parametrisation, along the diagonal of the pion-kaon mass plane we 
estimate m^ = 110 ±20 MeV (see Fig.[TJ). At high values of mx the phase boundary is not 
renormalisation scale dependent. We can clearly see the scaling region of the trictitical 
point as well. Unfortunately we were not able to locate it directly because the model 
parametrisation breaks down just before arriving at the TCP along the m^ = line, but 
based on the scaling near TCP we estimate it to be in the interval m K CP € {1700, 1850} 
MeV. We can directly compare our result on the location of the TCP with the recent 
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Figure 1. Phase boundary of the LcrM obtained using a one-loop parametrisation at T = 
and the formulae of large N c ChPT for continuation on the — m^-plane. 

lattice result of Ref. [ ^] by showing the phase boundary in the quark mass-plane. In 
Ref. [[121 the physical point is more close to the boundary than in our case and the value 
of the mass of the strange quark at the TCP is mj cp = 3m s (mf ys is the value of the 
strange quark mass at the physical point). In our case mj cp = 13 — 15 x m^ hys . In Ref. [ 
IT2"] the location of the TCP was estimated also by fitting with the corresponding scaling 
equation (m U) d ~ (mj cp — m s ) 5//2 ), but only points with m s < mf ys . The investigation 
performed in the LaM shows that the scaling region sets in very close to the m u = 
(m n = 0) axis. It would be interesting to know how the lattice estimate would change if 
points more closer to the m u = axis were available. 
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